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ABSTRACT 

o 
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£S| , We investigate the formation by accretion of massive primordial protostars in the 

range 10 to 300 Mq. The high accretion rate used in the models (M = 4.4 x 



C/3 



10" 3 M yr 1 ) causes the structure and evolution to differ significantly from those of 
both present-day protostars and primordial zero-age main sequence stars. After an ini- 
{Sj ' tial expansion of the radius (for M* < 12M©), the protostar undergoes an extended 

phase of contraction (up to M* ~ 6OM ). The stellar surface is not visible throughout 
most of the main accretion phase, since a photosphere is formed in the infalling en- 
velope. Also, significant nuclear burning does not take place until a protostellar mass 
(f-) , of about 80 Mq. As the interior luminosity approaches the Eddington luminosity, the 

protostellar radius rapidly expands, reaching a maximum around 100 Mq. Changes in 
the ionization of the surface layers induce a secondary phase of contraction, followed by 
a final swelling due to radiation pressure when the stellar mass reaches about 300 Mq. 
■ This expansion is likely to signal the end of the main accretion phase, thus setting an 

T^" 1, upper limit to the protostellar mass formed in these conditions. 

o 
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1. Introduction 

Numerical simulations are now beginning to constrain the properties of the sites where star 
formation first occurred in the post-recombination epoch (e.g., Abel, Bryan, & Norman 2000; 
Bromm, Coppi, & Larson 1999, 2001; Tsuribe 2001). These works indicate that the fragment mass 
scale of primordial gas clouds is relatively large (Mf ~ 10 3 Mq). On the other hand, the mass 
of protostellar cores formed inside these fragments is less than 1O _2 M0 at the beginning of the 
main accretion phase (Palla, Salpeter, & Stahler 1983; Omukai & Nishi 1998). Because of the large 
reservoir, the cores continue to grow in mass owing to accretion of ambient matter. Eventually, 
they become ordinary stars that shine by nuclear burning. The fundamental questions are then: 
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What fraction of the cloud mass can be incorporated into a star? Is there an upper limit to the 
mass of the primordial stars forming under such circumstances? 

In the context of present-day star formation, it is generally believed that stellar activity (ra- 
diation force, the expansion of an HII region, stellar winds, etc.) of sufficiently massive protostars 
halts accretion, thereby determining the final stellar mass (e.g., Larson & Starrfield 1971; Nakano, 
Hasegawa, &; Norman 1995). Of these processes, radiation force on dust grains appears to be the 
most efficient, making spherical accretion above a mass of ~ 15 M Q increasingly difficult for solar- 
type metallicities (e.g., Wolfire & Cassinelli 1987), although non-spherical accretion may enhance 
the maximum stellar mass (Nakano 1989). 

The main differences between primordial and Pop I accreting protostars are the following. First, 
the opacity of the infalling gas is drastically reduced due to the lack of dust grains. Second, the mass 
accretion rate, M, which is related to the sound speed c s in the protostellar cloud by M = c^/G 
(e.g., Stahler, Shu & Taam 1980), is substantially higher, because of the high temperature of the 
primordial clouds, typically T g ~ 1000 — 1500 K over a wide density range (Palla et al. 1983). 
These two properties imply a reduced effect of radiation pressure of the protostellar photons and a 
higher momentum of the inflow. Thus, reversal of accretion by radiation force alone seems to be less 
efficient, and a higher maximum stellar mass is expected (e.g., Larson & Starrfield 1971). Although 
qualitatively correct, this prediction has not been substantiated by numerical calculations covering 
a large mass interval. We still do not know the sequence of events that lead to the termination of 
the infall at the onset of protostellar surface activity. It is the purpose of this Letter to illustrate 
the results of such calculations for accreting protostars in the mass range 10 M & to ~ 300 M . 

2. Numerical Approach 

The evolution of low- and intermediate-mass primordial protostars has been followed by Stahler, 
Palla, & Salpeter (1986, hereafter SPS). In that study, an initially small core of 0.01 M Q grows 
rapidly by accretion at a rate of M = 4.4 x 10~ 3 M Q yr _1 , corresponding to the ambient gas 
temperature. As shown by Omukai & Nishi (1998), this is consistent with the physical conditions 
in cloud fragments at the time of core formation. The calculation of SPS was terminated at a mass 

=10.5 M Q , when nuclear burning had not yet started and the core was well below the mass 
scale of the parent cloud. At that point, the protostar was undergoing an internal readjustment 
due to strong heat transfer, resulting in the appearance of a luminosity wave at the surface and 
radius swelling. 

The basic strategy and equations used here are the same as in SPS. In this scheme, the proto- 
stellar evolution is treated as a sequence of steady state accretion flows onto a growing hydrostatic 
core (see Fig.l of SPS for a definition of the various regions). The core is assumed to be in hydro- 
static equilibrium and the ordinary stellar structure equations are applied. If the gas in front of 
the accretion shock is optically thin, no envelope model is constructed and the boundary condition 
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eq. (4a) of SPS is applied at the core surface. In case the preshock gas is optically thick, we 
integrate the equations inside the radiative precursor from the photosphere to the core surface (see 
eqs. 7a-d of SPS). Outside the photosphere, we assume a free- falling flow and evaluate at each 
time the deceleration exerted by radiation pressure. We adopted the same mass accretion rate of 
M = 4.4 x 10~ 3 M yr _1 as in SPS. The main difference with that work is in the use of updated 
opacity tables: for T < 7000 K, we take the results of Lenzuni, Chernoff, & Salpeter (1991, compo- 
sition of X = 72, Y = 28), and the OPAL opacity at higher temperatures (Iglesias & Rogers 1996, 
with slightly different composition X = 70, Y = 30). 

3. Results 

The evolution of primordial protostars in the range O.OIMq to 1O.5M0 has been computed by 
SPS. According to their results, early in the evolution, the accreted material settles adiabatically. 
Because of the high value of M, the core is surrounded by an optically thick radiative precursor. 
As the protostar grows, the luminosity increases until the Kelvin-Helmholtz (KH) time scale be- 
comes shorter than the accretion time. The contraction causes the outward propagation of the 
accumulated heat as a luminosity wave, which reaches the surface at a core mass of M* ~ 8M & . 
As a result, the radius starts increasing, and the radiative precursor disappears. Owing to the 
complexity of the internal readjustment caused by the luminosity wave, we decided to follow again 
this phase and chose the M* ~ 8 M G protostar of SPS as the initial model. 

3.1. Overall evolution 

The evolution of the core (R*) and photospheric (R p h) radii is presented in Figure 1. Initially, 
the expansion of the core continues from M* =8 to 11.5 M Q , reaching a maximum value of 220 R Q . 
Then, due to the very short KH time, the protostar starts contracting, and the radiative precursor 
appears again at a mass of M* = 12.4 Mq. Therefore, except for the short interval between 
M* ~ 8 and 12 Mq, the optically thick precursor persists throughout the main accretion phase. 
Unlike the core radius, the photospheric radius steadily expands as the photospheric luminosity 
(L p h) increases (see Fig. 2). The spatial extent of the precursor roughly corresponds to that of 
the ionized region because the photospheric temperature, T c g ~ 6000 K, is close to the ionization 
temperature. However, the present case is different from an ordinary HII region in that here the 
gas is optically thick to electron scattering and H~ bound-free absorption, in addition to Lyman 
continuum. 

The evolution of the interior (L*) and photospheric luminosities divided by the core mass is 
shown in Figure 2. The photospheric luminosity is the sum of L* and L acc , the accretion luminosity. 
As the luminosity-to-mass ratio L/M* grows, the effect of radiation force increases relative to that of 
gravity. Also, the accretion luminosity L acc ~ GM*M / R* increases as a result of core contraction. 
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Because of this effect and the high accretion rate, L p ^ approaches the Eddington limit (^Edd) at 
M* ~ 60 Mq. At this point, the core radius has shrunk to R* ~ 12 Rq, a value which is still higher 
than that of a zero-age main sequence star (2.8 Rq for a 50 M star, Marigo et al. 2001). 

Despite the large increase of the interior luminosity, accretion is not halted by the effect of 
radiation pressure. What happens instead is a second phase of expansion of the surface regions of 
the core starting at M* ~ 60 Mq (see Fig. 1) and a concurrent reduction of the accretion luminosity. 
The reason for the rapid increase of the radius is the following. While the inner part of the core 
continues contracting (see the curve of i?*(90%) in Fig. 1), the interior luminosity L* rises roughly 
as Ml and becomes close to the Eddington limit at M* ~ 80 Mq. At this time, hydrogen burning 
begins and core contraction stops. Owing to the high interior luminosity, strong radiation forces 
act on the thin surface layers, where the opacity sharply peaks due to the ionization of hydrogen 
and helium. This causes an acceleration of the expansion of the core from ~ 80 Mq to ~ 100 Mq. 
Then, the ratio L*/M* levels off at a value close to, but still smaller than, the Eddington limit (see 
Fig. 2) and the radius stops increasing. 

From then on, the surface layers of more massive protostars are subject to a delicate balance 
between the radiation force and the external pressure of the accreting envelope, whereas the interior 
relaxes toward the structure appropriate to a main sequence star (see the evolution of i?*(90%) 
in Fig. 1). The interior luminosity is always close to the Eddington limit and then increases 
approximately as M*, although the ratio L^/M* keeps increasing slowly. After further shrinking of 
the core (between 100 and 260 Mq) , the photospheric luminosity eventually reaches the Eddington 
limit. The physical conditions in the surface layers are thus similar to those that occurred at the 
time of the first expansion, but now the course of events is much more dramatic. The interior 
luminosity L* continues to increase and eventually equals the Eddington limit at ~ 300 M Q . This 
causes the huge expansion of the outermost layers of the protostar shown in Fig. 1, while the 
interior does not contract anymore. The core runaway signals the end of the accretion phase. 

3.2. Evolution of the Core Interior 

At the beginning of the calculations, the temperature is highest off-center (around q = M/M* ~ 
0.1) and the innermost part of the core remains adiabatic. In addition to this thermally inactive 
part, the core consists of a radiative region and a superadiabatic surface layer. The newly ac- 
creted material is first incorporated into the superadiabatic layer, where the material is heated by 
radiation. In the radiative relaxation layer, the material loses entropy and contracts. 

Because of the off-center temperature maximum, deuterium burning starts off-center at about 
12 Mq. The deuterium luminosity reaches a maximum between 12 and 16 Mq, but the contribution 
to the total luminosity is always negligible (<10%). Therefore, unlike Population I and II protostars, 
in primordial objects deuterium burning does not play an important role. The central inert part 
persists until the stellar mass reaches ~ 40Mq , when radiative heating from outside becomes strong 
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enough and the off-center temperature maximum disappears. We can see the corresponding sudden 
rise (drop) of the central temperature (density) in Figure 3b. Thereafter, the central temperature 
continues to rise gradually due to KH contraction. 

When the core mass reaches ~ 80 M@, the central temperature exceeds 10 8 K, needed to ignite 
the 3a reactions. Then, enough carbon is built up to drive the CN cycle and significant hydrogen 
burning ensues (see the sharp rise of L nuc in Figure 2). The onset and spread of central convection 
is shown in Figure 3a. Although carbon continues to be synthesized at the center, the growth of 
the convective core and mixing limits its abundance to values ~ 10~ 9 (see Figure 3b). 

Before the onset of the CN-cycle, the luminosity is generated mostly by contraction, with a 
minor contribution of H-burning via the pp-chain. This situation differs from the standard ZAMS 
models studied by several groups who find that the CN-cycle dominates the energy budget above 
~ 30Mq (e.g. Ezer & Cameron 1971). Since the high entropy of the accreted matter prevents 
contraction, the critical mass for the CN-cycle is much higher in the protostellar case. 

The radiation flow is absorbed by the high opacity in the superadiabatic layer. Because 
of the accumulated radiation energy in these regions, radiation pressure exceeds gas pressure at 

~ 60 Mq. The radiation dominated region quickly expands as the core mass increases. On the 
other hand, in the central convective core, the entropy increases in time due to nuclear burning. As 
a result, the central temperature remains almost constant, whereas the density decreases gradually 
(see Figure 3b). The ratio of radiation to total (gas+radiation) pressure keeps increasing, and 
eventually the central part becomes dominated by radiation at = 145 Mq. This region expands 
rapidly and merges with the outer radiation dominated region at 180 Mq. In more massive objects, 
the internal structure is fully determined by radiation pressure. 

4. Discussion 

As a result of the high value of the mass accretion rate characteristic of primordial gas clouds, 
the structural properties of massive protostars have several unique and unexpected features. The 
evolutionary tracks of the photosphere and core surface in the HR diagram are displayed in Figure 
4. The locus of the population III ZAMS stars computed from static models is also shown (Marigo 
et al. 2001; Bromm et al. 2001). The effective temperature of the protostellar core is computed 
from the core radius and preshock luminosity. The striking feature of the photospheric track is 
the almost constant value of the effective temperature. The photosphere forms in the accreting 
envelope, and persists almost throughout the main accretion phase. Because of the sensitivity of 
the H~ bound-free opacity, the temperature is locked at the value of T e $ ~ 6000 K (see eq. 23b of 
SPS). Observationally, primordial protostars of vastly different luminosities should have the same 
optical colors. On the other hand, the core surface will remain unobservable from outside because 
of the presence of the optically thick precursor. 

From Fig. 4 we see that the track of the core surface cuts through the ZAMS track of the core 
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surface at a mass ~ 40 M . However, the curve does not join the ZAMS smoothly, but follows a 
more elaborate path due to the behavior of the mass-radius relation. Therefore, realistic ZAMS 
models of Pop III stars should start with initial conditions that differ from those usually adopted 
on the basis of thermal equilibrium. 

Also, the nuclear properties derived from the protostellar phase are quite distinctive. Deu- 
terium burning and the p-p chain reactions have little effect on the evolution. Hydrogen burning 
due to the CN-cycle begins only at ~ 80 M . At this mass, the interior luminosity becomes close 
to the Eddington limit, so that the surface regions are subject to strong radiation forces. These 
conditions are favorable to the onset of activity in the form of mass loss. It is clear that continuum 
radiation driven winds can occur in Pop III protostars. According to the mass-radius relation of 
Fig. 1, there are two critical stages when the conditions become favorable for the onset of a wind: 
the rapid expansion of the surface layers at M* ~ 80 and 300 M . However, it is only at the highest 
mass that we expect a major episode of mass loss, since earlier on the interior luminosity always 
remains slightly below the critical value. Even in this case, however, we anticipate a short lived 
phase of wind activity related to the small extent of the high opacity surface layer. 

Our result of a upper mass limit of ~ 300 M is sensitive to the adopted value of the mass 
accretion rate and to its time evolution. In fact, M may not be constant and decrease in time from 
an even larger initial value (e.g., Omukai & Nishi 1998). Quantitatively, it is difficult to predict 
the exact values of M™ ax as a function of M without detailed models. However, we can expect 
that the evolution obtained here will not change qualitatively. In case of larger values of M, the 
higher luminosity will cause an earlier stripping of the surface layers by radiation pressure. Also, 
nuclear burning will be postponed at higher masses because of the shorter contraction time scale. 
On the contrary, for lower accretion rates, the evolution would resemble that of Pop I protostars 
with the central core smoothly becoming a ZAMS star during the main accretion phase, but with 
a fundamental difference. Because of the lack of grains, the radiative force is not capable to reverse 
the infall until a very large value of the protostellar mass. Thus, it is evident that there must be a 
critical value of M for the onset of the runaway increase of the protostellar radius that determines 
the maximum accumulated mass. The existence of this critical behavior will be explored in a 
forthcoming paper. 

This work is supported in part by Research Fellowships of the Japan Society for the Promotion 
of Science for Young Scientists, grant 6819. K.O. wishes to thank the director of the Osservatorio 
di Arcetri for the hospitality during an extended visit in 2001. 
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Fig. 1. — Mass-radius relation for massive primordial protostars. The evolution of the core (solid 
line) and photospheric (dashed line) radii is shown as a function of the protostellar mass. The 
circle on the curve of i?* marks the initial model of our calculations. The i?* and R p h curves at 
smaller masses come from SPS. The dotted lines represent the radii of mass shells containing 90% 
and 99%, respectively, of the total stellar mass. 
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Fig. 2. — Evolution of the luminosity to mass ratio as a function of mass. The interior and photo- 
spheric luminosities are displayed by the solid and dotted lines, respectively. The solid horizontal 
line represents the Eddington luminosity for electron scattering. The nuclear luminosity due to 
H-burning is shown by the dashed line. 
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Fig. 3. — The internal structure of the core as a function of mass, (a) Evolution of the nuclear energy 
and convection as a function of the relative mass q = M/M*. Regions where the energy generation 
by D- and H-burning exceeds 10% of the average energy generation rate L^/M* are shown as shaded 
regions The extent of the convective core is illustrated by the short-dashed area. The curve labelled 
P T = P g is the locus where the radiaton pressure equals gas pressure. Radiation pressure dominates 
to the right of this curve, (b) Evolution of central density, temperature, and mass fraction of 
CN-elements. Each variable is normalized to 0.1 g cm~ 3 , 10 7 K, and 10 -9 , respectively. 
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Fig. 4. — HR diagram for primordial protostars. The evolution of the photosphere and core surface 
is shown by the thick solid and dashed lines, respectively. For comparison, we also show the locus 
of the metal- free ZAMS stars (Marigo et al. 2001 for M* < 100 M ; Bromm, Kudritzki, & Loeb 
2001 for higher mass). The numbers on both tracks label the value of the core mass (in solar units). 
The filled circles on the dashed line have the same meaning. 



